Measuring Polynomial Invariants of Multi-Party Quantum States 
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We present networks for directly estimating the polynomial invariants of multi-party quantum 
states under local transformations. The structure of these networks is closely related to the structure 
of the invariants themselves and this lends a physical interpretation to these otherwise abstract 
mathematical quantities. Specifically, our networks estimate the invariants under local unitary (LU) 
transformations and under stochastic local operations and classical communication (SLOCC). Our 
networks can estimate the LU invariants for multi-party states, where each party can have a Hilbert 
space of arbitrary dimension and the SLOCC invariants for multi-qubit states. We analyze the 
statistical efficiency of our networks compared to methods based on estimating the state coefficients 
and calculating the invariants. 



I. INTRODUCTION 

Entanglement is a key resource in quantum informa- 
tion and computation since it can be used to perform 
tasks such as teleportation, super-dense coding and key 
distribution. Therefore, it is important to find ways 
of classifying and quantifying the entanglement prop- 
erties of quantum states. Central to this is the idea 
that locally invariant quantities can be used to char- 
acterize entanglement. Invariants under Local Uni- 
tary (LU) and more general transformations, such as 
Stochastic Local Operations and Classical Communica- 
tion (SLOCC), have been extensively studied in this con- 

text Bisaiiaaiiamiii. 

However, invariants are rather abstract mathematical 
objects and it is natural to ask whether any physical 
meaning can be given to them. One way of doing this 
is to investigate how these quantities might be measured 
given a number of copies of an unknown state. This 
could be done by simply measuring the coefficients of 
the state and then calculating the invariants. However, 
finding procedures to measure the invariants directly may 
be more efficient and also lends the invariants a physical 
interpretation as "collective observables" of the state. 

For bipartite pure states, the Schmidt coefficients are 
a complete set of LU invariants and optimal protocols 
for measuring them were given in |13| . Also, in [bH ] a 
method was given for estimating the polynomial SLOCC 
invariants of a general two-qubit state. 

In this paper we present networks for estimating two 
classes of polynomial invariants for multi-party states: 
the LU invariants for multi-party states with arbitrary 
local Hilbert space dimension and the SLOCC invariants 
for multi-qubit states. In both cases, the protocol works 
for both pure and mixed states. In particular, the struc- 



ture of the networks reflects the structure of the invari- 
ants in a very simple way. 

In Jill we review the construction of local invariants 
under LU transformations and in illlll the networks for 
measuring these invariants are presented. We then turn 
to invariants under SLOCC transformations, reviewing 
their construction in fllVl and presenting networks to mea- 
sure them in fJ2] In order to construct the networks for 
SLOCC invariants we make use of the Structural Physical 
Approximation (SPA) to non-physical maps introduced 
in 15]. The relevant details of this are presented in ^VII 
Finally, in flVIII we evaluate estimation protocols based 
on our networks by comparing them to simple techniques 
based on estimating the state coefficients. 



II. POLYNOMIAL INVARIANTS UNDER LU 
TRANSFORMATIONS 



A. Pure states 

Two n-party pure states , \t/j') € 
equivalent under LU transformations if 
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where Uj G U(dj) is a unitary operation acting on the 
Hilbert space of the jth party. States on the same orbit 
under this action have the same entanglement properties. 
Given a particular state, we might be interested in de- 
termining which orbit it belongs to. This can be done by 
establishing a canonical point on each orbit, such as the 
Schmidt form for bipartite states. However, canonical 
forms rapidly become more complicated as the number 
of parties is increased. Alternatively, we can construct 
polynomial functions of the state coefficients that are in- 
variant on each orbit. Theorems from invariant theory 
guarantee that a finite set of such polynomials is enough 
to distinguish the generic orbits under this action. We 
now review the construction of such a set. 
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FIG. 1: Diagrammatic representation of the quartic two-qubit 
LU invariant J, given in eq.ljlj. The first index of each term 
is represented by a circle and the second by a square. A line 
joins indices that are contracted with a 8. 



1. One party 

Consider the state \tp) — Ei=i a% I*) m a single party 
Hilbert space C d , where {\i}} is an orthonormal basis. 
The only independent invariant under unitary transfor- 
mations of this state is the norm (ip\ tp). This may be 
written as 



a" Si a 



(2) 



t-,3 



where S\ is the Kronecker delta. 8f is the U(d) invariant 
tensor and invariants for larger numbers of parties are 
formed by similar contractions of the state coefficients 
with their complex conjugates. 



2. Two qubits 

As an example, consider a two-qubit state l^) — 
El j=o aU Ni)- There is only one independent quadratic 
invariant, which is simply the norm of the state. How- 
ever, at quartic order we find the following invariant, 
which is functionally independent of the norm 



J 



= E< 



(3) 



For two qubits, we know that this is the only other in- 
dependent invariant because every state has a canon- 
ical Schmidt form \tp) = y/p\00) + ^/T^p\ll), with 
1/2 < p < 1 and J = 2(p 2 — p) + 1 determines p uniquely. 

Another useful way of representing the invariant is to 
define two permutations a, r on the set {1, 2} where a is 
the identity permutation and r(l) = 2, r(2) = 1. Then 



•W)=E ailil 



a J a, „ a, 

*<r(l)Jr(l) *<rC2)J-r(2) 



(4) 



This also suggests a diagrammatic way of representing 
the invariant (see fig. UJ. 



A multipartite pure state can be written in terms of 
an orthonormal basis as follows 
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i,j,k... 
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ijk... 



ijk. 



(5) 



A general polynomial function of the state coefficients 
and their complex conjugates can be written as 
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If the polynomial @) has equal numbers of a's and a*'s 
and all the indices of the a's are contracted using the 
invariant tensor 5 with those of the a*'s, each index be- 
ing contracted with an index corresponding to the same 
party then the polynomial is manifestly invariant under 
LU transformations. 

Such polynomials can be written in terms of permuta- 
tions on the indices. Let r be the degree of the polynomial 
in a (and hence also the degree in a*). Let <7,t, //... 
be permutations acting on the set {1,2, ... ,r} and let 
a = (c, r, /i, . . •)• Then the invariants can be written as: 



v iijiki.. i 2 j2k 2 
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(7) 



In fact, a can always be chosen to be the identity per- 
mutation by permuting the a terms in this expression. 
Additionally, each Jg can be associated with a diagram 
constructed in the same way as figH 

The invariants Jg are enough to completely distinguish 
the generic orbits under LU transformations. In fact, 
invariant theory guarantees that only a finite collection 
of them are needed to do this. However, except in a 
few simple cases, it is unknown which Jg invariants form 
minimal complete sets. 



B. Mixed states 

Two mixed states p, p' are equivalent under LU trans- 
formations if 

p' = Ux ® U 2 ® . . . ® U n pU\ ® Ul <g> . . . ® Z/t (8) 

The LU invariants for mixed states can be derived by 
rewriting the pure state invariants Q in terms of the 
density matrix p = {ip\ and noting that the resulting 
expressions are still invariant under LU transformations 
for general density matrices. This can be done by noting 
that terms such as Oi n ^"'Q.*j are elements of the den- 
sity matrix. A general density matrix may be written in 
terms of an orthonormal basis as 
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ij k... 
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(9) 
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FIG. 2: General construction of network to measure polyno- 
mial LU invariants. 



and the corresponding expression for an LU invariant is 



J_ _ _*ljlfcl— iafrkl... 
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III. MEASURING INVARIANTS UNDER LU 
TRANSFORMATIONS 

A. Network construction 

The general construction of the network used to mea- 
sure the LU invariants is shown in fig El It generalizes 
networks for estimating functionals of bipartite states 
given in [lil liH To measure an LU invariant of 

degree r in a (and also degree r in a*) we take r copies 
of the unknown state p. In addition, we take a single 
qubit in the state |0) and apply a Hadamard rotation H 
to transform the state to (|0) + |1)). In the next step, 
we apply a unitary operation U on the r copies of p con- 
trolled by the Hadamard rotated qubit. Finally we per- 
form a measurement on the single qubit in the {|0) , |1)} 
basis. The expectation value of this measurement will be 



(Z) = Re(Tr(t/p® r )) 



(11) 



When p — is a pure state then this is equivalent 

to 



(12) 



In order to determine networks for measuring the LU 
invariants, it only remains to show that there is a U such 
that the invariants can be expressed in the form (|llfl . 

To do this for pure states, we have to express poly- 
nomials of the form (JJJ in the form of (|12fl . Firstly, we 



note that a Ijfe " = (ijk. . .| ip), a 



ijk.. 



ijk. 



and 



to each permutation a in @ we associate a permutation 
matrix 
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FIG. 3: Network for measuring the 2-qubit quartic invariant. 



where P a acts on the Hilbert space of the same party for 
each of the r copies of the state \ip). Then to each a we 
associate the permutation matrix 



P« = P« 



Pt®Pu 



(14) 



where P a , P r , P„, . . . act on the Hilbert space of the same 
party as a,r,fx,... in Q on each of the r copies of the 
state. Then (J7J) can be written as 



j s = {i>f r p s W 



(15) 



Since Pg is unitary these invariants can be estimated with 
the network in fig[21by setting U = Pg to obtain the real 
part and U — iPg to obtain the imaginary part. For the 
specific example of the 2-qubit invariant © we have 



Ji 



(°-.t) 



A 2 B 2 * MM 



E) SWAP BlB2 



(16) 
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Note also that the physical construction of Pg is closely 
related to the diagram associated with Jg (compare figs, 
^and 13 for example). 

Finally, note that if p is a mixed state then applying 
the same procedure without modification will give the 
invariants of eq. ((TU1) . 

It has previously been noted 0] that all homogeneous 
polynomial LU invariants are determined by the expecta- 
tion values of two observables on r copies of a state. Here, 
we have given an explicit network for measuring these 
observables. Also, similar constructions can be made to 
estimate other polynomial functionals of quantum states 
[l7| and these can be modified to enable the estimation to 
proceed by LOCC [3, i.e. with no collective operations 
over the n-parties. A similar modification would enable 
the LU invariants to be estimated by LOCC, but this 
would affect the efficiency of the estimation discussed in 

mm 



IV. POLYNOMIAL INVARIANTS UNDER 
SLOCC 

When attempting to classify entanglement, it is of- 
ten useful to consider invariants under local transfor- 
mations that are more general than unitary transforma- 
tions. For this purpose, invariants under SLOCC have 
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been introduced^]. In 0we construct a network to 
measure the modulus squared of these invariants for the 
case where each party has a single qubit (i.e. the Hilbert 
space is (C 2 )®"). 



A. Pure states 

Two n-party pure states and are equivalent 
under SLOCC if it is possible to obtain \ip') with non- 
zero probability via a sequence of Local Operations and 
Classical Communication (LOCC) starting from a single 
copy of and vice- versa. In |20|. this criterion was 
shown to be equivalent to 



|^') = Mi ® M 2 ® . . . <8> M n \%l>) 



(17) 



where Mj € GL(dj) is an invertible linear transformation 
acting on the dj -dimensional Hilbert space of the jth 
party. 

In what follows, we find polynomial invariants for the 
special case where Mj £ SL(2), i.e. the transformation 
has unit determinant and each party has a single qubit. 
Networks to determine the modulus squared of these in- 
variants will be given in Note that it is not pos- 
sible to measure the SL(2) n invariants directly because 
they are not invariant under global phase transforma- 
tions \ip) — ► e 10 \ip), which have no physical significance. 
It is for this reason that we instead measure the modulus 
squared, which is invariant under these phase transfor- 
mations. 

Under general GL(2) n transformations, the polyno- 
mial SL(2) n invariants are still invariant up to a mul- 
tiplicative factor, which is just some power of the deter- 
minant of Mi ® M2 <8> • • • <8> M n . Thus, ratios of appropri- 
ate powers of these polynomials will be invariants under 
GL(2) n . 



1. Two qubits 

In order to illustrate the polynomial invariants under 
SL(2) n , first consider the case where n = 2. Two states 

= £j,*=i a 3 '* and W) = Ei,*=i «° fe \j k ) satisfy 
(O if 



a' = MyolMI 



(18) 



This means that det(a) = det(a') is an SL{2) x SL{2) 
invariant, since det(Mi) = det(M2) = 1. This may be 
written as 



(19) 



where the totally antisymmetric tensor is the SL(2) 
invariant tensor. For two qubit pure states, this is the 
only independent SL(2) x SL(2) invariant. 



2. General case 

The SL(2) n invariants can be constructed in a similar 
way to the LU invariants except the invariant tensor is 
now eij, and we contract a's with a's instead of a*'s. 
Thus, polynomials of the form 



Kg — 'Y2i\ e iii2 e jlj2 e kik 2 ■ ■ ■ e i 



-lir e j r -ljr e k,— 
Q ,j < T(l)J T (l)fc M (l)--- Q ,V(2)i T (2)fcf I (2)- 



(20) 



are manifestly invariant. Note that it is straightforward 
to generalize this construction to the case where each 
party has a <i-dimensional Hilbert space by contracting 
with the SL(d) n invariant tensor e^^...^ instead of e^. 
However, it is not yet clear how to measure these invari- 
ants because the effect of the higher rank e tensors cannot 
be physically implemented by linear transformations on 
states. 



B. Mixed states 

In general, two mixed states p, p' are equivalent un- 
der SLOCC if there exists two completely positive maps 
£1,82 which are implemcntablc via LOCC with non- 
zero probability of success such that p' = £\{p) and 
p = £2{p')- In order derive invariants using the expres- 
sions from the previous section, we will restrict to the 
case where p and p' are related by 

p' = Mi ® M 2 ® • . . O M nP M\ <g> M\ (g) . . . <g> Ml (21) 

with Mj € SL(2). The resulting expressions may not be 
invariant under more general SLOCC transformations, 
but are related to important quantities in entanglement 
theory as described in illV CI 

Unlike the LU invariants, it is not clear that H2(J|) can be 
written simply in terms of the coefficients of the density 
matrix p 
follows 



|?/>) ( , i/>|. However, \Kg\ can be written as 



\K,\ 2 = ^ 2 



2_,l e iii2 e jlj2 € k 1 k 2 



ljr e k r 



mim 2f nii! 2f piP2 



Vi) Ji-CDfenCl)- M2)Jt(2)*7.(2)- 
f"'<r(l)''T(l)I'/l(l)-/ J "»»(2)»T(S)P)l(2)- 

r™.(r)"T(r)P P W- 



(22) 

and these will also be SL(2) n invariants for mixed states. 



C. Examples of SL(2) n invariants 

The Kg invariants are especially interesting in entan- 
glement theory because many important entanglement 
measures can be easily calculated from them. For ex- 
ample, in the case of two-qubits, the concurrence |2lJ 
is defined as a simple function of the eigenvalues of pp, 
where 



P 



><7yP (Jy 



(23) 
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FIG. 4: Diagrammatic representation of the 3-tangle. The 
first index of each term is represented by a circle, the second 
by a square and the third by a triangle. A line joins indices 
that are contracted with an e. 



and T stands for transpose in the computational basis. 
These eigenvalues can be calculated from Tr((pp) m ) for 
m = 1, 2, 3, 4, which are simply the moduli squared of Kg 
invariants. In ^4j, networks were constructed to estimate 
these invariants for two qubits and we will generalize this 
construction to Kg invariants for larger numbers of par- 
ties. 

Another interesting example is the 3-tangle [22L l23l |. 
which is defined for pure states as the modulus of the 
following 3-qubit Kg invariant. 



T3 
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The 3-tangle gives information about the genuine 3-party 
entanglement between the qubits. 

Finally, note that the Kg invariants can be given sim- 
ilar diagrammatic representations to the Jg invariants. 
This is illustrated for the 3-tangle in figlH 



V. MEASURING SLOCC INVARIANTS 

The modulus squared of the SLOCC invariants can be 
measured using a network similar to fig. |2 except that 
the unknown states p must be preprocessed prior to the 
controlled-?/ operation. If Kg is of degree r in a then 
we will need r copies of p. The preprocessing stage will 
consist of collective unitary operations and completely 
positive maps that act on the entire Hilbert space of the 
r copies of p. The resulting state p', will yield the expec- 
tation value 



(Z) =Re(Tr(C/p')) 



(25) 



for the measurement at the end of the network. In this 
section, we describe the preprocessing operations and 
unitary operations U that enable the modulus squared 
of the SLOCC invariants to be written in this form. 

First, we apply the inverse of the permutation matrix 
associated with a to the r copies ofpto obtain Ptp® r Pg. 

The second, and final, part of the preprocessing stage 
is to apply a completely positive map A to the state. To 
describe A we first define the multi-party analogue of eq. 
(G3>. 



<7yp <Jy®(Jy< 



(26) 



Next, we define a map A that acts on a product of r states 
by applying the tilde operation to the even numbered 
states as follows 



A(pi ® p2 



Pr) — Pi ® P2 <8 P3 



' Pr (27) 



where each pj is an n-party state. 

Unfortunately, A cannot be physically implemented, 
since it is not a completely positive map. This can be 
dealt with by using the Structural Physical Approxima- 
tion (SPA) to A, which we will call A. A is the "closest" 
physical map to A. This is discussed in 3 VII but for now 
we construct the network as if A could be implemented 
perfectly. 

The final pre-processed state p 1 will be 



K{Plp^Pg) 



(28) 



Next, the controlled-C/ operation in our network must 
be chosen such that (Z) — \K |i when p' is used as the 
input. One can easily verify that the pairwise SWAP 
gate, defined by 



u\4>i)®\4>2) 
\<h)®\<h) 



<g> \4> r - 
•<8> \<j> r ) 



<g> \<j> r ) 
' \<f>r-l) 



(29) 



where \<f>j) is an n-party state fulfils this condition. 



VI. THE STRUCTURAL PHYSICAL 
APPROXIMATION 

The A operation encountered in the previous section 
is an example of a positive, but not completely positive 
map. These cannot be implemented exactly, but instead 
we can apply an approximation. 



A (p) = al + (3A(p) 



(30) 



where I is the identity operator and a, (3 are real positive 
constants chosen such that A is completely positive. If 
we fix a and (3 such that A is trace-preserving and f3 is 
maximized, then the results of |T^ imply that 



Hp) 



2*' 



i 



2 § nr + 1 2* 



23' 



1 



AGO (31) 



where n is the number of qubits in each copy of the state 
and r is the degree of the Kg for which we are estimating 
the modulus squared. 

On replacing A with A in our network the expecta- 
tion value of the Z measurement still allows the modulus 
squared of the Kg invariant to be determined via 



\Kg\ 2 = (2>" + l) (Z) 



Z)-T 



(32) 



However, the SPA does affect the accuracy to which the 
invariant is determined. This is discussed further in the 
next section. Additionally, in [l^, it is shown that this 
sort of SPA can be implemented by LOCC. Thus, the 
SLOCC invariants could also be estimated by LOCC, 
but the efficiency discussed in WII Bl would be affected. 
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VII. EVALUATION 



A. Statistical analysis of the network 



The main aim of the protocols presented in gTTTjand gV| 
is to provide a physical interpretation for the polynomial 
invariants. However, we have not yet addressed the ques- 
tion of how efficient these measurement protocols are. In 
this section, we compare the efficiency of our protocols 
to protocols based on simply measuring the state coef- 
ficients and calculating the invariants. We use unbiased 
estimators based on counting j^lSEIlBl- Also, we per- 
form the analysis in the limit where a large number of 
copies of the state have been measured, so that the vari- 
ances of the estimates are small and can be treated to 
first order in all subsequent calculations. We note that 
more sophisticated estimation procedures are also possi- 
ble |27| , but our purpose here is to compare the networks 
to methods that are easily accessible experimentally. 

Measuring the state coefficients would clearly be a 
more straightforward procedure to perform experimen- 
tally than using our network. Although more parameters 
have to be determined, this does not necessarily mean 
that it is a less efficient method for estimating the invari- 
ants than using our networks. There are several quite 
general reasons why this might be the case. 

Firstly, suppose that we are interested in measuring 
a complete set of polynomial LU invariants for some 
unknown state of n parties, where each party has a d- 
dimensional Hilbert space. In general, we do not know 
how many we would need to measure, but parameter 
counting arguments Q, 0, show that the number of 
local degrees of freedom is linear in n whereas the total 
number of degrees of freedom is exponential in n. Thus, 
for large n almost all the degrees of freedom are non- 
local. Even for moderately sized n, there are nearly as 
many invariants as there are state coefficients. In addi- 
tion, the invariants are typically highly non-linear func- 
tions of the state coefficients. For these reasons, we ex- 
pect that measuring a complete set of invariants directly 
will generally not be more efficient than measuring the 
state coefficients for large n. Similar considerations also 
apply to the SLOCC invariants. 

Despite these considerations, it may be the case that 
our networks are more efficient if we are only interested 
in measuring a small incomplete subset of the invariants. 
Also, they may be more efficient for estimating complete 
sets when n is small. For this reason, and for simplicity, 
we concentrate on estimating two qubit invariants in this 
section. 

There are also other reasons why our protocols may 
not be efficient. For example, our protocols only em- 
ploy a two-outcome measurement for each r copies of 
the state whereas estimating the state coefficients uses a 
two-outcome measurement on each copy. Also, for the 
Kg invariants, we will sec that using the SPA introduces 
a lot of noise into the measurement. Nonetheless, there 
are still some cases where using our networks is more 
efficient than estimating the state coefficients. 



For a particular setup in our network we make repeated 
measurements of an observable Z ', wit h exp ectation value 
F = Tr {Up'). Z is a random variable [31) with distribu- 
tion 



p(Z = +l) = Ul + F) 
p(Z = -l) = Ul-F) 



(33) 



If we define the event Z = +1 as a success and set 
p = P(Z = +1) then repeating the network N times 
is equivalent to performing N Bernoulli trials. The num- 
ber of successes N s is a random variable with a binomial 
distribution and its expectation value is (N s ) — Np = 
Tr (1 + F). In an actual experiment, the observed num- 
ber of successes N s can be used to compute an unbiased 
estimator for F, given by 



N s 



1 



with variance 



var 



1 

N 



(34) 



(35) 



We are interested in determining how many trials are 
needed in order for the estimate F to be reasonably ac- 
curate. Specifically, we would like to quantify how many 
trials are needed to make the variance of var(F) < e 
for some e > 0. In an experimental situation, we would 
not be able to calculate var(F) from our data, so we 
would have to estimate it using the sample variance, 
var(F). However, in the limit N — > oo we can use the 
fact that var(F) = 0(N~ r ) and var(var(F)) = 0(7V" 4 ), 
i.e. vhr(F) converges to the true variance much faster 
than F converges to F so var(F) var(F). Thus, in this 
limit we have that 



N > -(l-F 2 ) 
e 



(36) 



Recall that for the LU invariants, the real and imagi- 
nary parts of the invariant are estimated independently 
and that each use of the network requires r copies of the 
state, where r is the degree of the invariant in a. If we 
use the same number of samples for estimating both the 
real and imaginary parts then the total number of copies 
required is 



M>-(2-\J ff \ 
e 



(37) 



In some cases, we know a priori that the invariant is 
always real or always imaginary. If this is the case, then 
we can achieve the same accuracy with 



M > - (1 



(38) 



For the SLOCC invariants, each use of the network 
requires r copies of the state, where r is the degree of the 
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invariant in a. Also the estimate of the invariant must 
take into account the use of the SPA via H32[) . In this 
case, the total number of copies required is 



M > 



{2i nr + lf - (\K ff \ 2 + T 



(39) 



Notice that the 2 3nr term will dominate the term in the 
square bracket for large n and r. This is due to the noise 
introduced into the measurement by the SPA. 



B. Comparison to methods based on state 
estimation 

In order to evaluate our protocols, we compare them 
to methods based on estimating the density matrix of 
the state and then calculating the invariants. We do this 
by estimating each state coefficient using observations on 
single copies of the state. This is known as homodyne to- 
mography (see [23] for an overview and also 0EaiE3)- 
This is not the optimal way of reconstructing the state 
in general psj. but it will greatly simplify the analysis. 



1. Example: Two-qubit LU invariants 
A general two-qubit density matrix can be written as 
P=\ {li®h + a j a j ®h + 6,^2 ® &j 



(40) 



The two-qubit LU invariant @ can be written in terms 
of these coefficients as 



J = Tr04) 



(41) 



Each bj can be determined by simply performing a 
<jj measurement on Nj copies of Bob's half of the state. 
The probability distributions of the associated random 
variables are given by 



p(<7,- = +L) = i (1 + 6,0 
p(aj = -1) = § (1 - 6,0 



(42) 



Thus, each bj can be estimated in the same way as F 
in 1|34|) and we have that 



var(6j) 



(43) 



We can then construct an estimator for J given by 



J 



(44) 



which will be biased, but in the large Nj limit 



var( J) 




(45) 



to first order in var(6j). 

If we make the additional restriction that each ob- 
servable o j is sampled the same number of times (i.e. 



Nj = y) then we must take 



(1 



(46) 



for our estimate to have variance J$ e. 

One way to compare this to the result for our net- 
work is to take an average over all pure states. If wc 
assume that all pure states are equally likely, i.e. inte- 
grate (IVII Bl and gnj using Haar measure (for details 
see |29jL then we find that on average we will need 3/2 
times as many copies of the state if we use the coefficient 
estimation method. This is half of what one might ex- 
pect from parameter counting alone, since three times as 
many parameters are estimated in the state coefficient 
method. The factor of two is explained by the fact that 
each use of our network uses two copies of the state. 

However, it is possible to find parameter ranges in 
which the state coefficient method performs better than 
our networks. One such range is given by setting 61 = 



0, 



< 63 < 



This illustrates the fact that 



parameter counting does not always reflect the statistical 
efficiency of a given protocol. Any partial information we 
have available about the type of states being measured 
might change our judgement of which protocol is more 
efficient. 



2. Example: Two-qubit SLOCC invariants 

For the two-qubit SLOCC invariants we take the 
quadratic invariant l|19l) as an example. In terms of the 
decomposition (|40|) this can be written as 



\K\ 2 = 



1 



(47) 



If we estimate this by measuring all 15 of the state coeffi- 
cients an equal number of times then by a similar analysis 
to the LU case we find that we need at least 



6,(1-6?)] 



R%) 



(48) 



copies of the state to get a variance ;$ e. 

Taking averages, one finds that fewer copies are needed 
in the state coefficient protocol by a factor «5x 10 3 de- 
spite the fact that many more parameters have to be 
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estimated in this protocol than when using our network. 
This is largely due to the factor 2 12 that appears in (|39|l . 
which arises from the noise introduced by the SPA. This 
suggests that other estimation and detection protocols 
based on the SPA [LJ, [l6j may be less efficient than pa- 
rameter counting arguments would imply. In fact, there 
are no states for which our network performs better than 
the coefficient estimation method. Even in the best pos- 
sible case for our network, the state coefficient method 
requires fewer states by about 3 orders of magnitude. 

VIII. CONCLUSIONS 

We have presented networks for measuring the polyno- 
mial invariants of quantum states under LU and SLOCC 



transformations. The structure of these networks is 
closely related to the structure of the invariants them- 
selves and thus gives the invariants a physical interpreta- 
tion. Comparison of these networks with methods based 
on estimating the state coefficients indicate that the net- 
works are of limited practical use for estimating complete 
sets of invariants. Indeed, our results suggest that any es- 
timation procedure that employs the SPA is statistically 
inefficient even when the number of parties is small. 

We know that no procedure for estimating invariants 
directly can outperform protocols based on estimating 
the state coefficients as the number of parties is increased. 
For small numbers of parties it seems that there can be 
some increase in efficiency, but the optimal protocol is 
not known in general. 
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